UNITUE-THEP-01 /OO 
hep-ph/0001261 



A Comparison between Relativistic and Semi— Relativist ic 
Treatment in the Diquark— Quark Model^ 



M. Oettel 1 and R. Alkofer 2 

Institute for Theoretical Physics, Tubingen University 
Auf der Morgenstelle 14, D- 72076 Tubingen, Germany 



Abstract 

In the diquark-quark model of the nucleon including scalar and axialvector di- 
quarks we compare solutions of the ladder Bethe-Salpeter equation in the instan- 
taneous Salpeter approximation and in the fully covariant (i.e. four-dimensional) 
treatment. We obtain that the binding energy is severely underestimated in the 
Salpeter approximation. For the electromagnetic form factors of the nucleon we 
find that in both approaches the overall shapes of the respective form factors are 
reasonably similar up to Q 2 ~ 0.4 GeV 2 . However, the magnetic moments differ 
substantially as well as results for the pion-nucleon and the axial coupling of the 
nucleon. 
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Bound states of two particles in quantum field theory can be described by solutions of 
integral equations such as the Bethe-Salpeter (BS) equation |J. It has the generic form 

*(p,P)=5 x (pi)5 2 (p 2 ) J^K(p,p',P)^(p',P), (1) 

where Si and S2 denote the fully dressed propagators of particle 1 and 2, and K is the full, 
2P irreducible interaction kernel between the in- and outgoing particles. The wave function 
(or Bethe-Salpeter amplitude) ^> is not an observable quantity but it can be employed to 
compute e.g. form factors as matrix elements of operators 0. 

Little success has been made in going beyond the ladder approximation (see, however, 
e.g. ||) where the propagators Si and S 2 are approximated by the corresponding bare 
propagators and the kernel K just contains the interaction to lowest order which might be 
the exchange of a third particle with bare propagator S3: 

*(p, p) = s b r e ( P i) sr e {P2) j ^^(p,j/, p)w, n (2) 

Even in this approximation an analytical solution to the BS equation has been found only 
for the Wick-Cutkosky model with a bound state of zero mass (for an introduction into 
this model see e.g. Q). Moreover, practical attempts to solve the BS equation in various 
models have resorted to further approximations, especially three-dimensional reductions 
(for a critical comparison of different three-dimensional reductions, see ||). The Salpeter 
approximation in which retardation effects in S3 are neglected reduces the BS equation to 
a Klein-Gordon or Dirac equation with a hermitian Hamiltonian. It is often claimed to be 
covariant though the neglection of the time-component of the four-momenta p and p' seems 
to prohibit this phrasing. It is true, however, that the solutions obey relativistic translation 
invariance. 

In this letter we choose a diquark-quark model of the nucleon for a comparison between 
the Salpeter approximation and the fully relativistic BS solution. This has two reasons: If 
one assumes separable diquark correlations, the three-body Faddeev equations reduce to a 
ladder BS equation which sums up the quark exchange between the spectator quark and 
the diquark quasiparticle || [j], ||, ||, [10], [TT| . Thus the truncation of the BS kernel to the 
one-particle exchange is not completely arbitrary but directly linked to the separability as- 
sumption. Secondly, a comparison between relativistic and semi-relativistic solutions in the 
BS equations should not stop at the level of the solutions to the equations but involve observ- 
able quantities. There are several attempts to calculate observables within this framework, 
such as the electromagnetic form factors. This has been done in the Salpeter approximation 
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]i~2 , p~3f] and using the full solution Jl4 , 15[| . In the following we will compare the electromag- 
netic form factors of the nucleon obtained in the Salpeter approximation to the ones using 
the full BS equation employing model parameters as in refs. [12|, [L3| . 



The diquark-quark model 

The underlying approximations of the diquark-quark model are, first, the neglection of 
any 3-particle irreducible interactions between the quarks in the nucleon, and, secondly, the 
separability of correlations in the two-quark channel. The first assumption allows to derive 
a relativistic Faddeev equation for the 6-point quark Green's function, and the second one 
reduces it to an effective quark-diquark BS equation. 

The full derivation of the quark-diquark BS equation within the NJL model can be 
found in refs. |7], |9|. In the following we will sketch the most important steps. As usual in 
the treatment of the BS equation we will work in Euclidean space. Scalar and axialvector 
diquarks are introduced via the assumption of a separable 4-point quark function: 

P is the total momentum of the incoming and the outgoing quark-quark pair, p and q are 
the relative momenta between the quarks in the two channels. Xapip) and x^p{p) are vertex 
functions of quarks with a scalar and an axialvector diquark, respectively. They belong 
to a 3-representation in color space and are either flavor antisymmetric (scalar diquark) or 
flavor symmetric (axialvector diquark). For their Dirac structure we will retain the dominant 
contribution only. Thus we introduce one scalar function P(p) which depends only on the 
relative momentum p between the quarks, 

Xap(p) = g s (l 5 CUP(p), (3) 
X^(P) = 9a (YCUP(p). (4) 

C denotes hereby the charge conjugation matrix, and g a and g s are effective coupling con- 
stants obtained by normalization of the diquark states. The Pauli principle requires the 
relative momentum to be defined p = \{p a — Pp), where p a and pp are the quark momenta 
|p~5f . P(p) parametrizes the extension of the vertex in momentum space. To facilitate the 
comparison between our work and that of [|T^, we choose 



P(p) = exp(-4AV) • (5) 
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The propagators of scalar and axial vector diquark are the ones for a free spin-0 and spin-1 
particle, 

D (p) = ~ 2I 2 > ( 6 ) 

P + m sc 

D^( P ) = 2 \ ■ (7) 

p A + m Q 



ax 



£ is a gauge parameter introduced in |]16| , and in the following we will put £ = 1. The 
constituent quark propagator is simply the free fermion propagator 

Sip) = ■ (8) 

p z + m z 

The nucleon BS wave function can be described by an effective multi-spinor character- 
izing the scalar and axialvector correlations (see e.g. [|TT| ) , 



where u(P) is a positive-energy Dirac spinor with P being the total momentum of the 
bound state, p is the relative momentum between quark and diquark, respectively. The 
vertex function is defined by amputating the legs off the wave function, 

* n ° dm. do) 



The BS equation is a system of equations for wave and vertex function that takes the 
compact form 

P ■G(p,p',P) L)(P',P) = (11) 



(2tt) 4 ' 

where the object G(p, p', P) involves the propagators of quark and diquark and the interaction 
kernel that describes the quark exchange between quark and diquark, 



G(p,p',P) = (27rr5(p-p')S- 1 (p a )^ D Q 1 {D ^ yi y Pd 



if xS T (q)x -V3 X »'S T (q)x 

2 \-V3 X s T {q)r -x»'s T ( q )r 



The partitioning of the momentum between quark and diquark introduces a parameter 77 
with p q = rjP + p and pd = (1 — rf)P — p. Therefore the momentum of the exchanged quark 
is q = —p — p' + (1 — 2t])P. The relative momentum of the quarks at the diquark vertex x is 
P2 = p+p'/2— (1— 3r])P/2 and the one at the conjugated vertex x ispi = p/2+p'— (1— 3rj)P/2. 
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Relativistic translation invariance requires that if $(p, P; rji) is a solution of the BS equation 
then also $(p + (772 — T)i)P, P\ V2) is one. 

/ n \ 

In this 
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The BS equation (|Tl| ) is solved in the rest frame of the bound state, P - 
frame the Salpeter approximation amounts to neglecting the fourth component of all vectors 
appearing in the interaction, i.e of the vectors q, pi and p2- The immediate consequence is 
that the vertex function will depend only on the relative three-momentum, P) = 



Numerical Solutions 

In the following we use the complete partial wave decomposition of wave and vertex 
function for octet baryons given in | |1 1| | , 



/ 



$ 5 (p,P) 
$ 4 (p,P) 



P 



\ 



V 



1A 2 
ip(a-p)A 3 + (cr x p)(o-p)A 5 
- P A 4 + 1((tx p)A 6 



(13) 
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The unknown scalar functions Si and Ai are functions of p 2 = p^p^ and of the angular 
variable z = P ■ p which denotes the cosine of the angle between p^ and the 4-axis. As 
explained above, the functional dependence collapses in the Salpeter approximation from 
the two variables {p 2 ,z) to the single variable p 2 (l — z 2 ). We expand the scalar functions in 
terms of Chebyshev polynomials of the first kind in the variable z, 



SAA^p 2 ^) = J2^ n S?[A7](p 2 )T n (z), 



(14) 



n=0 



and derive a system of coupled integral equations for the Chebyshev momenta Sf [Af] that 
we solve iteratively as described in [11 , |I~5| . This expansion is very close to the hyperspherical 
expansion that has been shown to work extraordinarily well in the massive Wick-Cutkosky 
model |T7| , |l~8f and in quenched QED [Tj|. This finding has been corroborated by |TTJ in the 
diquark-quark model for point like diquarks. 

As stated already we adopt the parameters of refs. |12|, [TJ). In the case of the scalar- 
axialvector diquark model this especially includes identical diquark normalizations leading 
to equal diquark-quark-quark couplings in the scalar and axialvector channel, g s = g a . The 
BS equation ([□]) can be written as an eigenvalue problem for g s . This coupling is adjusted 
to yield the physical nucleon mass M=0.939 GeV for given values of the quark and diquark 
mass as well as the diquark width A. In table |] we have listed the two sets of parameters, 
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Set 




m q 
[GeVJ 


m sc 
[GeVJ 


[GeVJ 


A 
[rmj 


9s 




I 


Salpeter 
full 


0.35 


0.65 




0.18 


20.0 
16.5 


(20.0) 


II 


Salpeter 
full 


0.35 


0.65 


0.65 


0.24 


12.3 
9.6 


(11.5) 



Table 1: The two parameter sets of the model. 



Chebyshev Moments of S, Chebyshev Moments of S, 



Salpeter approximation Full calculation 




1 2 3 4 5 ' 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

p [GeV] p [GeV] 



Figure 1: The Chebyshev expansion of the dominat scalar function S\(p, z) in the Salpeter 
approximation and the full calculation. 

scalar diquark only and scalar-axialvector diquark model, with their corresponding eigen- 
values obtained by us in the full calculation and the Salpeter approximation. The values 
in parentheses are the ones from refs. JT2"| , [TB|. Please note that due to a different flavor 
normalization these values had to be multiplied by to be directly comparable to ours. 

Although we could reproduce the eigenvalue for the model case with scalar diquark 
only, this is not the case for Set II. We observe that the calculations of [TB| involved only 



4 instead of 6 axialvector components of namely the projected ones onto zero orbital 
angular momentum and the corresponding lower components. Still one would expect a higher 
eigenvalue in the reduced system. The more striking observation is the amplification of the 
eigenvalue by about 20... 25% in the Salpeter approximation although the binding energy is 
small, being only 6% of the sum of the constituent masses. This is in contrast to results 



obtained in the massive Wick-Cutkosky model |L9| where the Salpeter approximation leads 
to a reduction of the eigenvalue. This may be attributed to the exchange of a boson instead 
of a fermion as in the present study. 

The substantial difference between the two approaches is also reflected in the vertex 
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Figure 2: Impulse approximation diagrams contributing to the electromagnetic current. 

function solutions themselves. Fig. [I] shows the Chebyshev moments of the dominant scalar 
function Si for both methods using the parameters of Set I. Only the even momenta are 
given, since the odd ones are zero in the Salpeter approximation (but are present, of course, 
in the full calculation). Two things are manifest: The Salpeter amplitudes have a much 
broader spatial extent than the full amplitudes. Secondly, the expansion in Chebyshev 
polynomials that relies on an approximate 0(4) symmetry converges much more rapidly for 
the full solution but is hardly convincing in the Salpeter approximation. Again, since the 
scalar functions depend in the Salpeter approximation on just one variable, p 2 (l — z 2 ), our 
expansion is a cumbersome way of visualizing the solution but makes clear that the Salpeter 
approximation strongly violates the 0(4)-symmetry. 

Electromagnetic Form Factors 

The Sachs form factors Ge and Gm can be extracted from the solutions of the BS equations 
using the relations 

G E =^- 2 Tr(J»)P», G M = ^Tr(J^ } (15) 

<^> = fPu f pt- A s> T {p f ,Pf)j»m,Pi)- (i6) 



q(p,P n )=MA+ (17) 

P = Pn 



, (2tt) 4 J (2tt) 4 

with the definitions P = (Pi + -P/)/2 and 7^ = 7 M — P^ff . To this end, one has to normalize 
the wave and the vertex function through the condition 

r d A P r d A P ' - 
~J (2^7 W) " 

which uses the definition of G given in eq. (|12|) and employs the positive-energy projector 
A + . This normalization integral is again performed easily in the rest frame of the bound 
state, additionally the double integral over the interaction kernel drops out in the Salpeter 
approximation since the kernel is independent of P in the rest frame. 

The current operator J M in eq. (|1(J) contains all possible couplings of the photon to 
the kernel G of the BS equation. The two simplest contributions make up the impulse 
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approximation pictorially shown in Fig. |2|. Since the impulse approximation does not mix 
scalar and axialvector amplitudes there are altogether four diagrams to compute: quark and 
diquark coupling and for each of the two the matrix element between scalar and axialvector 
amplitudes of the nucleon. Flavor algebra [13| yields the following current matrix elements 
for proton and neutron, 

2 1 

(Jproton) = ~^(Jq,s) + ~^(^dq,s) + (^dq,a) ' \^°) 

im 1111 

{Jneutron) = ~~^(Jq,s) + g(^d?,s) + gWg/i) _ gWd?,a)' (^) 

The diagrams of the impulse approximation separately conserve the current in the Salpeter 
approximation. The proof depends on the behaviour of the current under time reversal and 
parity transformation and can be found in jl2| and more explicitly in ||20|| . Furthermore 



there is the peculiar identity at zero photon momentum transfer, 

(^ sW W 2 = 0)) = (^, s[a] (g 2 = 0)), (20) 

which guarantees that proton and neutron have their correct charge (i.e. G P E (Q 2 = 0) = 1, 
G^(Q 2 =0)=0). This is not so in the full calculation, to ensure current conservation and the 
physical charges one has to take into account all diagrams of the photon coupling to the 
interaction part of G. These diagrams are shown in fig. |3|, the proof of this assertion and the 
construction of the "seagull" diagrams (photon coupling to the diquark-quark vertices) can 



be found in see also |]2T| for a general discussion of the current operator in three-body 
theory. 

In fig. |] the electric form factors of proton and neutron are displayed, using parameter 
set II. The first observation is that in the Salpeter approximation we could not obtain 
convergence with the expansion in Chebyshev polynomials beyond Q 2 ~ 0.4 GeV 2 . We 
computed the form factors in the Breit frame where Q is real but Zi = pi ■ Pi and Zf — pf Pf 
have imaginary parts and their absolute values may exceed one, except for the case of no 



momentum transfer (TTJ. So this expansion that works in the rest frame, and it does barely 
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Figure 3: Exchange diagrams contributing to the electromagnetic current. 
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Figure 4: Nucleoli electric form factors in the Salpeter approximation and in the full cal- 
culation. The curves for the Salpeter approximation have been obtained by including the 
Chebyshev moments of the vertex function up to order n max . 



so for the Salpeter approximation, will not generally work in a moving frame. On the 
other hand, the decrease of the higher Chebyshev moments is so drastic for the full four- 
dimensional solution that form factor calculations converge easily up to several GeV 2 [IB, ^2 . 
However, as can be seen from fig. |] the Salpeter approximation badly fails above 0.5 GeV 2 
thereby revealing its semi-relativistic nature. 

The second finding concerns the electromagnetic radii. The Salpeter approximation tends 
to underestimate the proton charge radius and to overestimate the absolute value of the 
neutron charge radius, see the first two rows of table The axialvectors tend to suppress 
the neutron electric form factor much more in the full calculation than in the Salpeter 
approximation. 

Turning to the magnetic moments, the contributions of the various diagrams are tab- 
ulated in table for the proton. Following [13| we ascribed to the axialvector diquark a 
rather large anomalous magnetic moment of k—1.6 which was needed by the cited author 



Set 


^sc,q 


f^sc,dq 




f^ax,dq 




SUM 


I 


Salpeter 


1.57 


0.01 








1.58 (1.58) 




full 


2.04 


0.01 






0.33 


2.38 


II 


Salpeter 


0.92 (1.08) 


0.00 





1.37 (1.7) 




2.29 (2.78) 




full 


1.09 


0.00 





0.98 


0.38 


2.45 



Table 2: Contributions to the proton magnetic moment from the single diagrams. The first 
index refers to the nucleon amplitudes involved and the second to the photon coupling to 



either quark or diquark. In parentheses the values of refs. [12, 13] are given. 
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Set 




Dia- 






Hp 


ypjm 




{r n )m 


9tvNN 


T-itNN 


9a 








grams 


[fm] 


[fm 2 ] 




[fm] 




[fm] 




[fm] 




[fm] 


I 


Sal. 


imp. 


0.89 


-0.28 


1.58 


1.04 


-0.77 


1.05 


8.96 


1.04 


0.93 


0.93 




£ 11 

mil 


imp. 


1 r\r\ 

1.00 


-0.21 


2.05 


1.09 


1 AO 

-1.02 


1.09 


11 T 1 

11.71 


1 r\ r 7 

1.07 


1.16 


0.99 






all 


0.99 


-0.23 


2.38 


1.06 


-1.65 


1.00 


15.34 


1.04 


1.46 


1.02 


II 


Sal. 


imp. 


0.88 


-0.06 


2.29 


0.85 


-0.98 


0.94 


6.03 


1.15 


0.53 


1.08 




full 


imp. 


1.01 


-0.04 


2.07 


0.99 


-1.02 


1.06 


6.95 


1.21 


0.65 


1.14 






all 


1.01 


-0.04 


2.45 


0.99 


-1.26 


1.05 


9.36 


1.15 


0.82 


1.10 



Table 3: Some static quantities of the nucleon. The numbers for the Salpeter calculation 
and the ones in the first line for the full calculation are obtained with impulse approximation 
diagrams only. The second line for the full calculation includes the couplings to the quark- 
diquark interaction kernel. For g a and g n NN, only diagrams with quark couplings have been 
considered. 

to fit the proton magnetic moment. As already observed earlier, we could reproduce the 
magnetic moment for Set I, however, for Set II, the values differ and especially the coupling 
to the axialvector diquark is much weaker in our Salpeter calculation. Rather more inter- 
esting is the comparison in the full calculation between Set I and Set II: The axialvector 
diquark improves the magnetic moment only marginally! The Salpeter approximation tends 
to overestimate this contribution quite drastically. 

Finally we want to mention that the Salpeter approximation underestimates the pion- 
nucleon coupling g n NN an d the axial coupling g^ quite sizeably, see table ^. 

Conclusions 

In this letter we have presented results for a covariant diquark-quark model. The ladder 
BS equation for the nucleon has been solved in a fully covariant way and in the instantaneous 
Salpeter approximation. As for the model with scalar diquarks only we have verified the 



results of ref. [J2J whereas there are discrepancies if the axialvector diquark is included. Part 
of these differences are due to the fact that in ref. |H| not all (ground state) axialvector 
components have been taken into account. Additionally, we take our result as an indication 
that the calculations presented in ref. [|H| might suffer from some minor error. 

The main purpose of this letter is the comparison of observables calculated in the Salpeter 
approximation to the ones obtained in the fully four-dimensional scheme. The first very sur- 
prising observation is the overestimation of the BS eigenvalue in the Salpeter approximation. 
Phrased otherwise, for a given coupling constant the binding energy would be much too small 
in the Salpeter approximation. We have also demonstrated that the Salpeter approxima- 
tion violates badly the approximate 0(4) symmetry of the BS equation. This has drastic 
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consequences for the resulting nucleon electromagnetic form factors if the photon virtuality 
exceeds 0.4 GeV 2 . Whereas different nucleon radii differ only mildly in these two approaches 
one sees very clearly that the results for the magnetic moments, the pion-nucleon coupling 
and the axial coupling are underestimated in the Salpeter approximation. 
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